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Abstract 

Based on the rainbow approximation of Dyson- Schwinger equation and the assumption that 
the inverse dressed quark propagator at finite chemical potential is analytic in the neighborhood 
of [i = 0, a new method for obtaining the dressed quark propagator at finite chemical potential 
/x from the one at zero chemical potential is developed. Using this method the dressed quark 
propagator at finite chemical potential can be obtained directly from the one at zero chemical 
potential without the necessity of numerically solving the corresponding coupled integral equations 
by iteration methods. A comparison with previous results is given. 
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QCD at finite density is of fundamental importance, both on purely theoretical and 
phenomenological grounds. Although considerable progress has been achieved in the lattice 
treatment of finite density QCD [1-3] recently analytical studies with effective theories are 
useful alternatives and often shed more light on nonperturbative phenomena than the lattice 
approach. Due to the fact that the quark propagator at finite chemical potential plays 
an essential role in the study of chiral symmetry restoration and quark deconfinement, it 
is interesting to give a general recipe to study the chemical potential dependence of the 
dressed-quark propagator at non-zero chemical potential in the framework of a suitable 
nonperturbative QCD model. 

Since the Dyson-Schwinger(DS) approaches [4-6] provides a nonperturbative framework 
which admits the simultaneous study of dynamical chiral symmetry breaking and confine- 
ment, it is expected to be well suited to explore the transition from hadronic matter to 
QGP[7]. It is the aim of this letter to study the chemical potential dependence of the 
dressed-quark propagator in the framework of DS approach, which provides a means of de- 
termining the behavior of the chiral and deconfinement order parameters. Up to this end let 
us briefly review the DS approach used in this letter. The DS approach is based on a cou- 
pled set of integral equations between the quark, gluon, and ghost propagators and vertex 
functions. They form a countably infinite set of coupled integral equations with the equa- 
tion for an n-point Schwinger function depending on (n+1) and higher point functions. In 
order to handle this system it is necessary to make certain simplifications and truncations. 
One approach, which is commonly referred to as the "rainbow" approximation, employs 
the bare quark-gluon vertex and solves the Dyson-Schwinger equation(DSE) for the dressed 
quark propagator, with a given dressed gluon propagator g^D^^p) as input. For example, 
under the "rainbow" approximation, the quark self-energy E[/x](p) at the non-zero chemical 
potential /i can be written as 

4 r d 4 q 

SM(P) = 3 J j^u9s D ^(P ~ Qhv,G\M](Qhw (!) 
The quark propagator Q[fi](p) and the quark self-energy E[//](p) are related by 

£~V](p) = i7 • i>- ^74 + (p) 

4 f d A q 

= «7 • p - /174 + - J j^g 2 s D^(p - q)l£[n](q)lu- (2) 

Here we want to stress that Eq.(2), which employs the bare quark-gluon vertex and solves 
the DSE for the dressed quark propagator G\p](p) at finite chemical potential with a given 
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chemical potential independent gluon propagator g 2 D^(p) as input, is our starting point 
for studying the dressed quark propagator at finite chemical potential in the framework of 
DS approach. 

If /i is set to be zero, G[fj](p) goes into the dressed quark propagator in the zero chemical 
potential G{p) = Q\\i = 0](p), which reads 

G-\ V ) = i 1 ■ P A(p 2 ) + B(p 2 ) =ry- P + y (0I*#V(P - Qh?G(<lh», ( 3 ) 

with the self energy functions A(p 2 ) and B(p 2 ) being determined by the rainbow DSE in the 
chiral limit: 

for convenience, we have used a model ansatz g^D^p) = 5^ u g 2 D(p 2 ) for the gluon propa- 
gator, which is often referred to as the so-called "Feynman-like" gauge propagator [5]. It is 
readily seen that the B(p 2 ) in Eq.(4) has two qualitatively distinct solutions. The "Nambu- 
Goldstone" solution, for which B(p 2 ) 7^ 0, describes a phase, in which: 1) chiral symmetry 
is dynamically broken, because one has a nonzero quark mass function; and 2) the dressed 
quarks are confined, because the propagator described by these functions does not have a 
Lehmann representation. The other solution, the "Wigner" one, B(p 2 ) = 0, describes a 
phase, in which chiral symmetry is not broken and the dressed-quarks are not confined[4,5]. 
In the "Wigner" phase, the Dyson-Schwinger equation(4) reduces to: 

where A'{p 2 ) denotes the dressed quark vector self energy function in the "Wigner" phase. 

In order to motivate our new method for calculating the dressed quark propagator at finite 
chemical potential, let us first recall the methods adopted in the previous literatures[8-ll]. 
From Lorentz structure analysis, the most general form for the dressed quark propagator at 
finite chemical potential reads 

£~VKp) = *7 • PMp) + il4C(p) + B(p) - /i747 • PD(P), (6) 

where p = (p,p4 + ifJ,). Substituting Eq.(6) into Eq.(2) and by means of suitable projection 
procedure(multiplying by appropriate gamma matrices and then taking the trace), the four 



independent scalar function A(p), B(p), C(p), and V(p) are found to satisfy a coupled set of 
Dyson-Schwinger equations. 

In principle, for a given model gluon propagator, one can numerically solve these coupled 
DSEs. However, the numerical solution for an arbitrary gluon propagator turns out to 
be rather difficult in practice(In Ref.[10], the last term on the right hand of Eq.(6) was 
dropped). In order to avoid the difficulty of numerically solving four independent coupled 
integral equations, we have proposed a new method in Ref.[12] for calculating the chemical 
potential dependence of the dressed quark propagator. For the sake of completeness we 
outline the main points of this method here. Let us assume that the full inverse quark 
propagator at finite chemical potential is analytic in p, at least for small p. Under this 
assumption, one can expand (? _1 [/t](p) in powers of p as follows[12] 



dp, 



P + --- + 



i d n g- l [p] 



[i=0 



ni 



dp T 



p n + 



fi=0 



g- 1 + r (1 V + • • • + -r (n V + 



/i 



with rW(p, 0) 



rW(p,o) = 



d n Q-Mip) 



dp n 



n>\. 



(7) 



(8) 



fi=0 



It should be noted that Eq.(7) is only valid within the radius of convergence of p expansion. 

Working in the framework of rainbow approximation of DSE, one can prove the following 
by induction[12] 



r(")(p o) = dnG ~ 1(p) n>l 

Combining Eqs.(7) and (9) we have 

G-Mlp) = G- 1 (p) + r^(p,o)p + ... + ^ n \ P ,o)p n + 



(9) 



„_ lt , dG-\p) 
= G {p) + . , P + - 

= G-\p, Pi + ip) = G- 1 ( 1 



nl 

| 1 ^(p) | 
n\ d{—ipi) n 
J_ ^(p) 

+ n\ d{ Pi y m + 

= i~f-pA(p 2 ) +B{p 2 ). 



(10) 



Therefore the dressed quark propagator at finite p can be obtained from that at zero p by 
the simple substitution p 4 — * p 4 + ip. From this result one sees immediately that there are 
only two independent Lorentz structures(instead of four, which is determined from general 
analysis) in the dressed quark propagator at finite chemical potential. 



Substituting Eq.(lO) into Eq.(2) and then multiplying by appropriate gamma matrices 
and taking the trace in Eq.(2), we have the following: 



' d 4 q 2 _ A(q 2 )p ■ q 

(2tt) 4 ^ [P Q) g 2 A 2 (g2) + B 2 (q 2 ) 



B(q 



,~;2\ 



In principle, for a given model gluon propagator, one can numerically solve this set of coupled 
DSEs using iteration methods just as one solves the corresponding set of coupled DSEs at 
fj, = 0(Eq.(4)). However, in practice this is much more complicated than solving Eq.(4). Due 
to the presence of the \x in the fourth component of the momentum variable, the number 
of arguments of the dressed quark propagator and the independent integration dimension 
in Eq.(ll) is larger than those of Eq.(4) and this fact makes it difficult to find a stationary 
solutions by means of iteration method. In order to avoid this difficulty, here we shall adopt 
another more efficient method instead of the iteration method to solve Eq.(ll). This is the 
motivation of our present work. 

Applying the differential operation on both sides of Eq.(3), we obtain 



274 + 



d 4 q d 



g 2 s D(p-q) J u G{q)j u . 



(12) 



(2tt) 4 d(p 4 ) 

Similarly, applying the differential operation on both sides of Eq.(12) successively 
(n-l)(n > 2) times, we obtain 



d n G~\p) 
d(pA) n 



d 4 g d n 
(27r) 4 d(p,r 



g 2 s D{p-q) 7„G(qhv 



(13) 



Based on Eq.(lO) and Eqs. (12-13), we have main results in the present work 



8G^ 1 (n) 

g-Mip) = *7 • PA(P 2 ) + B(p 2 ) = G-\p) + n/ ^f V 4 



d(p 4 ) 







d 4 q 



1 d n 
n\ d(p 4 ) n 



g 2 s D{p-q) (in) n + 



g 2 s D(p-q)\ (^) + ...+ 
luG{q)-i u 



G-\p) -^4 + 



d 4 q 

W) 4 



g 2 D(p -q)- g 2 D(p - q) j u G{q)j v 



= i"f ■ p + 



d 4 q 
W) 



(14) 
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where we have made use of Eq.(3) in the last step in deriving Eq.(14). If \i is set to be 
zero, Eq.(14) reduces to Eq.(3). This is just what one would expect in advance. At this 
point a possible place which might lead to misunderstanding should be clarified. In our 
formulation, the ^-independent dressed gluon propagator g 2 D^ u (p) is taken as input and 
the dressed quark propagator at finite /i is obtained by solving the rainbow DSE(Eq.(2)). 
Therefore the explicit presence of // in g 2 D{p — q) does not imply that the gluon propagator 
becomes //-dependent. The //-dependent "gluon propagator" D(p — q) is only a quantity 
resulting from our mathematical trick(Eqs. (12-14)) and should not be interpreted as the 
actual // dependent gluon propagator. 

It should be noted that Eq.(14) has two major consequences. First, by means of Eq.(14) 
the dressed quark propagator at finite chemical potential can be obtained directly from the 
one at zero chemical potential without the necessity of numerically solving the correspond- 
ing integral equations by iteration methods(see Eqs. (15-16) below). This feature facilitates 
numerical calculations considerably. Second, Eq.(14) clearly shows that the whole nontriv- 
ial //-dependence of the dressed quark propagator is determined by the //-dependence of 
D(p — 9) (the 27 • p term only gives the "trivial" //-dependence). This means that, in a 
complete treatment of the dressed quark propagator at finite chemical potential, the depen- 
dence of the dressed gluon propagator on frequency cannot be neglected. Therefore, in all 
models adopting the instantaneous approximation (such as the NJL model) the nontrivial 
//-dependence of the dressed quark propagator will be lost (retardation effect are important 
here) . 

Multiplying by appropriate gamma matrices and then taking the trace in Eq.(14), we 
have the following: 

8 f d 4 q 2ri /~ A(q 2 )p-q 



[A(P } " 1]P = 3 J J2^ 9MP - q W(f) + B 2 {q 2 ) ' 
16 f d 4 q 2 B(q 2 ) 

Substituting B(p 2 ) = into Eq.(15), we have the "Wigner" solution for A'(p 2 ) 



The above result simplifies the calculation of A(p 2 ), B(p 2 ) and A'(p 2 ) greatly. Once the 
form of the model gluon propagator g 2 D(p — q) is given, one can determine A(p 2 ), B{p 2 ) 
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and A'(p 2 ) by solving the rainbow DSE (Eqs.(4) and (5)). Because the analytic expression for 
g 2 D(p — q) is available, one can determine g 2 s D{p — q) by analytic continuation, and thereby 
obtain A(p 2 ), B{p 2 ) and A'{p 2 ) using Eqs. (15,16) without the necessity of numerically solving 
the coupled integral equations of A(p 2 ) and B(p 2 ) by means of iteration method. From the 
above result one can easily obtain the dressed quark propagator at finite chemical potential 
from that at zero chemical potential. Here we want to stress that Eqs. (15) and (16) only 
holds under the "rainbow" approximation of DSE and within the radius of convergence of 
the fj, expansion. In the case of real QCD, it should be noted that both the dressed quark- 
gluon vertex and the dressed gluon propagator are chemical potential dependent. In this 
case, Eqs. (15) and (16) would fail. 

Just as was shown in Eqs. (15) and (16), the task of calculating the chemical potential 
dependence of the dressed quark propagator in the"Nambu-Goldstone" and the"Wigner" 
phase is reduced to the calculation of the three scalar functions A(p 2 ), B(p 2 ) and A'(p 2 ). In 
order to get the numerical solution of the above three scalar functions, one often use model 
forms for gluon two-point function as input in Eqs. (4) and (5). As a typical example, we 
choose the following model gluon propagator; 

g 2 D(q 2 ) = ATT 2 d—-——— with d = —. (17) 
y w ; q 4 + A 27 K J 

The dressed-gluon propagator in Eq.(17) simulates the infrared enhancement and confine- 
ment and it leads via the QCD gap equation(Eq.(4)) to an infrared enhancement of the 
light quark mass function. These modification are intimately related to the confinement and 
dynamical chiral symmetry breaking [4, 5]. The model parameters x an d A are adjusted to 
reproduce the weak decay constant in the chiral limit f n = 87 MeV. The forms of g 2 D(q 2 ) 
have been used in Ref. [13] and it has been shown that with these values a satisfactory 
description of all low energy chiral observables can be achieved. Substituting Eq.(17) into 
Eqs. (4) and (5), one can numerically solve the three scalar functions A(p 2 ), B(p 2 ) and A'(p 2 ). 
The calculated values of the above three scalar functions for model gluon propagator (17) 
with three sets of different parameters is plotted in Figs. (1,2). 

Based on the calculated values of A{p 2 ), B(p 2 ), A'{p 2 ) and Eqs. (15,16), it is not difficult 
to obtain the corresponding "Nambu" and "Wigner" solutions at finite chemical potential, 
i.e., A(p 2 ), B(p 2 ), A'(p 2 ). In order to check whether the numerical solutions A(p 2 ), B(p 2 ) 
obtained from Eq.(15) indeed satisfy the original coupled integral equations(Eq.(ll)), one 
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may directly substitute them into Eq.(ll) and we find that Eqs.(ll) are indeed satisfied. 
This result can be regarded as a self-consistent check of our numerical calculations in the 
present work. It is now apparent that our method has the merits that one can easily find 
the numerical solutions of A(p 2 ) and B(p 2 ) without the necessity of numerically solving Eq. 
(11) by iteration methods. 
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Fig.l. A(s) and A'(s) for gluon propagator (17) with three sets of different parameters . 
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Fig.2. B(s) for gluon propagator (17) with three sets of different parameters . 

With these two "phases" characterized by qualitatively different momentum-dependent 
quark propagator in the case of non-zero chemical potential, the DSE can be used to ex- 
plore chiral symmetry restoration and phase transition between the "Wigner" and "Nambu- 
Goldstone" phase. 

To explore the possibility of a phase transition one must consider the relative stability of 
the confined and deconfined phase, which is measured by the \i dependent vacuum pressure 
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difference (or "bag constant"), which is equivalent to calculating the difference between the 
tree- level auxiliary- field effective action [14] evaluated with the "Wigner" solution character- 
ized by B(p 2 ) = 0, and the "Nambu-Goldstone" solution characterised By B(p 2 ) ^ 0[8]: 



12 



/ 



dp 4 

W) 1 



In 



p 2 A 2 (p 2 ) + B 2 (p 2 ' 



p 2 A' 2 (p 2 ) 



+ Re 



p 2 A(p 2 ) 



p 2 A 2 (p 2 ) + B 2 (p 2 ) A'(p 2 ) 



(18) 



Substituting /i = into Eq.(18), we have the bag constant at zero chemical potential 
13(0). By means of numerical studies, the bag constant 13(0) for three different parameter 
sets of the model gluon propagator(Eq.(17)) are obtained and we list them in Tabled. 



Table. I. The bag constant for the model gluon propagator 





A [GeV 4 ] 


X [GeV] 


B^(0)[GeV] 


io- 1 


1.77 


0.122 


io- 2 


1.33 


0.126 


10~ 4 


0.95 


0.130 



The numerical results for the ratios of B(fi)/B(0) are plotted in Fig. 3. The scale is cal- 
culated to be B(0) = (0.122 GeV ~ 0.130 GeV) 4 (see Tabled), which can be compared 
with the value (0.145 GeV) 4 commonly used in bag-like models of hadron[15]. B(fi) > 
indicates the stability of the confined phase (Nambu-Goldstone) and hence the phase bound- 
ary is specified by £>(/i = /jl c ) = 0. B(fx) is positive when the "Nambu-Goldstone" phase is 
dynamically favored; i.e., has the highest pressure and become negative when the pressure of 
"Wigner" phase become larger. In Fig.3, we see that B(fi)/B(0) decrease with increasing fx 
within a certain range of the small chemical potential (It should be noted that our numerical 
results is only valid for small values of chemical potential. This is because that our model 
gluon propagator has no explicit fi -dependence while the actual gluon propagator should 
be \i dependent due to quark vacuum polarizations. As such it may be inadequate at large 
value of n). This result is qualitatively the same as that given in Refs.[7,10]. In Refs.[7,10], 
B(n)/B(0) decreases with increasing chemical potential up to /i c . 
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Fig.3. Ratio B(n)/B(0) as a function of fi. 

Now let us turn to the study of the measure of the dynamical chiral symmetry breaking 
in the case of non-zero chemical potential. In order to get a reasonable result for the mixed 
quark-gluon condensate and vacuum susceptibilities in an effective quark-quark interaction 
model, the authors in Refs. [16,17] defined the "effective" two-quark condensate as the dif- 
ference between the "exact" quark propagator(quark propagator in the "Nambu-Goldstone" 
phase, in which chiral symmetry is dynamically broken and the dressed quarks are confined) 
and the "perturbative" quark propagator (quark propagator in the " Wigner" phase, in which 
chiral symmetry is not dynamically broken and the dressed quarks are not confined). It can 
be written as (in the chiral limit and at zero chemical potential): 



<0|g#> M=0 = -tr DC {^[/i = 0] - QWfr = 0]} . 



(19) 



It should be noted that Eq.(19) is only valid in an effective quark-quark interaction 
model(more details can be found in Ref.[18]). 

Here we extend the above concept to get a measure of dynamical chiral symmetry breaking 
in the case of finite /i and obtain the "effective" two-quark condensate with the finite \i as: 



<0|w|0)„ = -tr DC {RegWfr] - ReQ^}} = -12 J ^-Re 



X2\ 



p 2 A 2 (p 2 ) + B 2 (p 2 ) 



.(20) 



Substituting /i — into Eq.(20), we have the usual "effective" two-quark condensate in 
the chiral limit. The calculated ratio (0|gg|0) M /(0|gg|0)^ =0 is plotted in Figs. (4,5). From 
Figs. (4,5) we can see that the behavior of (0|gg|0)^/(0|gg|0) M=0 depends strongly on the 
choice of parameters in the model gluon propagator. For the first set of parameters 
(0\qq\0) n/ (0\qq\0) n=o decreases monotonically as /x increases, while for the second and third 
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sets of parameters (0|gg|0) M /(0|gg|0) M= o first increases with increasing //, and after reach- 
ing a crest, decreases with increasing ji. This results is quite different from that given by 
Refs. [7,10,19-20]. In Refs.[7,10], one claims that (0|gg|0) M /(0|gg|0) M=0 increases with increas- 
ing \i up to /i c as a consequence of the necessary momentum-dependence of the dressed-quark 
self-energy. 




n(GeV) 



Fig.4. The ratio (01^^10)^/(01^^10)^=0 a s a function of /i. 
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Fig.5. The ratio (0|gg|0) M /(0|gg|0) M=0 as a function of /i. 



It should be noted that this strong dependence on model parameters is not unexpected 
because the whole nontrivial //-dependence of dressed quark propagator Q~ 1 [fj](p) is deter- 
mined by the //-dependence of D(p — g)(Eq.(14)). In practice this can be used as guide 
for choice of model parameters. Based on the Hellmann-Feynman theorem and Gell-Mann, 
Oakes, Renner relations(PCAC), one finds: 

(6|??|6>, 



(6|gg|6) 



= 1 



p=0 



+ 0(p 2 ). 



(21) 
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The leading-order linear baryon density(p) dependence is determined by the nucleon sigma 
term in a model-independent way(for details see Ref.[21]), i.e., <tat = {N\m q (uu + dd)\N). 
If we regard Eq.(21) as a criterion for studying the chiral order parameter at finite density, 
then one will find the second and third sets of parameters of the model gluon propagator (17) 
is not suitable. 

To summarize: in the present work, we provide a general recipe to calculate the chemical 
potential dependence of the dressed quark propagator under the rainbow approximation to 
the DS equation. This approach has the advantage that the dressed quark propagator at 
finite chemical potential can be obtained directly from the one at zero chemical potential 
without the necessity of numerically solving the corresponding coupled integral equations by 
iteration methods. This feature facilitates numerical calculations considerably. From this the 
"effective" quark condensate at finite chemical potential and the phase transition between 
the "Wigner" and "Nambu-Goldstone" phase is analyzed. It is found that the behavior 
of the "effective" quark condensate (6|gg|6) M /(0|gg|0) M=0 depends strongly on the choice of 
parameters of the model gluon propagator. This can serve as a guide for choosing model 
parameters in further study. Finally we want to stress that the approach adopted in this 
letter is general in principle and can be applied to the study of the color superconductivity 
in the framework of rainbow DS approach[22]. 
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